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Abstract. The kinematic expansion history of the universe is investigated by 
using the 307 supernovae type la from the Union Compilation set. Three simple 
model parameterizations for the deceleration parameter (constant, linear and abrupt 
transition) and two different models that are explicitly parametrized by the cosmic jerk 
parameter (constant and variable) are considered. Likelihood and Bayesian analyses 
are employed to find best fit parameters and compare models among themselves and 
with the flat ACDM model. Analytical expressions and estimates for the deceleration 
and cosmic jerk parameters today {qo and jo) and for the transition redshift (zt) 
between a past phase of cosmic deceleration to a current phase of acceleration are 
given. All models characterize an accelerated expansion for the universe today and 
largely indicate that it was decelerating in the past, having a transition redshift 
around 0.5. The cosmic jerk is not strongly constrained by the present supernovae 
data. For the most realistic kinematic models the la confidence limits imply the 
following ranges of values: go e [-0.96, -0.46], jo e [-3.2, -0.3] and G [0.36,0.84], 
which are compatible with the ACDM predictions, qo = —0.57 ± 0.04, jo — —1 and 
zt = 0.71 ± 0.08. We find that even very simple kinematic models are equally good to 
describe the data compared to the concordance ACDM model, and that the current 
observations are not powerful enough to discriminate among all of them. 

Keywords: supernova type la - standard candles, dark energy experiments, dark energy 
theory 



1. Introduction 



The extension of the Hubble diagram to larger distances by using observations from 
supernovae type la (SNIa) as standard candles, allowed the history of cosmic expansion 
to be probed with deeper detail. Independent measurements by various groups indicated 
that the current expansion is in fact speeding up and not slowing down, as believed for 
many decades [U El EJ H] . In other words, in virtue of some unknown mechanism, the 
expansion of the Universe underwent a "dynamic phase transition" whose main effect 
is to change the sign of the universal deceleration parameter q{z). 
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The physical explanation for such a transition is one of the greatest challenges for 
cosmology today. Inside the General Relativity paradigm, it requires the presence of a 
cosmological constant in the cosmic equations, or to postulate the existence of an exotic 
fluid with negative pressure (in addition to dark matter), usually called dark energy 
[3 El [3, El E] , or even a gravitationally-induced cold dark matter creation [101 El IE]- 
Another possibility is to change the theory describing the gravitational interaction as 
happens, for instance, in the framework of the so-called F(R) modified gravity theories 
[T3l [HI UHl HH [TTj. In both cases, the space parameter associated with the cosmic 
expansion is too degenerate, and, as such, it is not possible, based on the current 
data, to decide which mechanism or dark energy component is operating in the cosmic 
dynamics |T8]. 

Another very distinct firsthand approach to access the history of the cosmic 
expansion without the use of quantities coming from the dynamic description has also 
been proposed in the literature [I9]. This route is very interesting because it depends 
neither on the validity of any particular metric theory of gravity nor on the matter- 
energy content of the observed Universe. It is closely related to the weaker assumption 
that space-time is homogeneous and isotropic, so that the FRW metric is still valid, as 
are the kinematic equations for redshift /scale factor. Some call it cosmography [201 [22] 
or cosmokinetics [23], others use the term Friedmannless cosmology [2^ [25] , but, in 
what follows, we refer to it simply as a kinematic approach since it holds true regardless 
of the underlying cosmic dynamics [3l |26l [271 [2HI [291 jSO] ■ 

Few years ago, Elgar0y & Multamaki [25] investigated constraints on some 
kinematic models by employing a Bayesian marginal likelihood analysis based on the 
Gold Supernova sample data of Riess et al. [3] and the Supernova Legacy Survey (SNLS) 
of Astier et al. |4|. In their analysis of the fiat case, three different parameterizations 
for an accelerating q{z) model were examined, namely: constant, linear, and abrupt 
transition, respectively, Mq, Mi and M2 in their nomenclature. It was also argued that 
any expansion of the jerk parameter (the third order contribution in the expansion for 
kinematic luminosity distance in terms of the redshift z) could be seen as requiring more 
parameters, or cosmic fluids in the dynamic approach, than expanding the deceleration 
parameter. In addition, it was also observed that the flat ACDM and the Einstein-de 
Sitter model have constant jerk parameter jo = —1, and, therefore, it cannot be used 
for discriminating such cosmologies. Accordingly, their analysis was restricted to the 
deceleration parameter. 

In this work we go one step further by examining the case for the jerk parameter 
(constant and variable). One basic reason is that the bi-dimensional space parameter 
{%)jo) can naturally discriminate the fiat ACDM and Einstein-de Sitter models because 
they have different predictions for go- Actually, it is not necessary to expand i{z), 
thereby introducing many parameters in order to have a larger class of models, since this 
happens even for constant jerk. Probably, and more importantly, it is not clear a priori 
if the Bayesian evidence prefers the cosmic concordance model when a more general 
constant jerk parameter is considered. Potentially, a constant jerk parameter provide 
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us with the simplest approach to search for departures from the cosmic concordance 
model. For completeness, in our analysis we also consider all the flat models discussed 
by Elgar0y & Multamaki [25], however, we also examine the predictions of two different 
approaches including constant and variable jerk parameters. In addition, and differently 
from previous works, the present Bayesian analysis is based on the larger Union 
Compilation SNIa data recently published by Kowalski et al. [31] . 



2. Kinematic Models 



Assuming that the universe is homogeneous and isotropic above some scale, the 
Friedmann- Robertson- Walker (FRW) metric provides a good description of the geometry 
of the universe (c = 1) 

dr'^ 



ds' 



dr - a\t) 



+ r^dn 



(1) 



1 — kr"^ 

where k is the space curvature that we will assume hereafter to be null, since the current 
emerging consensus is that the universe is flat or very close to flat [18j. The function 
a{t) is the scale factor, which contains the complete history of the cosmic expansion and 
can be parameterized by the redshift, a = (1 + 2;)"^ (oq = 1). 

The rate of expansion and the acceleration are represented by the Hubble and 
deceleration parameters, respectively 
d 
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Similarly, the jerk parameter is defined as 
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The basic aim here is to examine some simple kinematic models for the cosmic expansion 
based on specific parameterizations for in ([3]) and a constant jerk parameter. 

The first and simplest model, Mq, is given by a constant deceleration parameter, 
qi^z) = go- The second model. Mi, is a linear expansion of the deceleration parameter 
q{z) = qo + qiz (first used by [3]). Model M2 (introduced by [19]) depicts two phases of 
constant deceleration parameter, separated by an abrupt transition redshift, q{z) = gg 
for z < Zt and q{z) = qi for z > Zt- The fourth model, M3, is a constant jerk 
parametrization, j{z) = Jq (examined for the first time by [3J). As one may show, 
the cosmic jerk is related with the deceleration parameter by the differential equation 

dq 



J 



q + 2q^ + {l + z] 



dz 



(5) 



It should be stressed that kinematic models with constant jerk parameter are also very 
attractive because the historically important dynamical EdS cosmology and, with more 
generality, the fiat ACDM scenario are particular cases for which j{z) = jo = —1. 
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On the other hand, it is widely known that if one wishes to describe the recent 
cosmic expansion, then the current values of the parameters given by ([2]), ([3]) and (jlj) 
lead to the late time cosmic expansion [32l [20] 

a{t) = 1 + H^{t - to) - \q^Hl{t - t^f - ^Jo^^o - ^o)' + 0[{t - to)% (6) 

from which the luminosity distance can be expanded, yielding an extended version of 
the Hubble law [33] 

+ 0{z% (7) 

where the highest order term depends on fourth order and higher derivatives of the scale 
factor. A shortcoming of this series is that for supernovae at ^ > 1 the 0{z^) terms 
may be large in principle, i.e., there is a convergence problem [3l]. Nevertheless, most 
of the known supernovae are at z < 1 and a truncation of ([7]) always can be seen as a 
polynomial fit. So we define our last model, M4, as the third order truncation of ([7]), 
which also has as free parameters go and jo- However, differently from M3, this model 
has variable jerk parameter (see IA.191) . 

In the Appendix, one may find the basic analytical expressions for all models 
investigated in the present work. 



dAz) 
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3. Sample and Statistical Analysis 

In the statistical analysis below we consider the most complete supernovae data set 
currently available, namely, the Union supernova sample as compiled by Kowalski et 
al. [31]. The Union SNIa compilation is a new data set at low and intermediate 
nearby-Hubble-flow redshifts whose analysis procedures permit to work with several 
heterogeneous supernova samples. It includes 13 independent sets, and, after selection 
cuts, the robust compilation obtained is composed by 307 SNIa events distributed over 
the redshift interval 0.015 < z < 1.62. 

In figure 1 we show the behaviour and dispersion of the data in the reduced Hubble- 
Sandage diagram. The curves correspond to the models considered in this work, as 
indicated in the legend. For the sake of comparison, we also show the present cosmic 
concordance model. 

The distance moduli /i is defined as the difference between the apparent (m) and 
absolute (M) magnitudes, so that the observed and theoretical values are respectively 

^^obs,^ = mobs,i - M , (8) 

and 

fith{zi) = mth{zi) - M = hlogiQ diizup) + /iq , (9) 
where /io = 25 — 5 log^^g -^o, and di is the luminosity distance 

dL{z,p) = {l + z) Tt^, (10) 

io H{u) 




Figure 1. Kinematic model predictions for SNe la data. Residual magnitude versus 
redshift is displayed for 307 SNe type la from SCP union compilation. The predictions 
of five kinematic models {Mq — M4) are displayed relative to an eternally coasting 
model, for which H{z) = Ho{l + z) and q{z) = 0. For comparison, the cosmic 
concordance model has also been included. Note that for z > 1.5, the model with 
j is a free parameter constant (M3) is the closest one to ACDM. 



which carries the model parameter dependencies represented by p. 
The hkehhood analysis is based on the calculation of 



[fJ'obs,i — fJ'th{Zi)y 
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We analytically marginalize over the nuisance parameter /io 



/ + CO 
exp 
-00 



to obtain 



where 



X =a + — 

c Vivr 



SNIa 
SNIa 



[Slogio di{Zi,p) - fiobs,'. 

51ogio di{zi,p) - fiobs,i 



(11) 
(12) 

(13) 

(14) 

(15) 
(16) 



Bayesian Analysis and Constraints on Kinematic Models from Union SNIa 6 



SNIa ^« 



The nuisance parameter value that minimizes f|T2|) is yUo = ^/c The expression 
X^{Pih/c) = a — {b'^/c) is sometimes used instead of f[T^ to perform the hkehhood 
analysis. Both are equivalent if the prior for /io is flat, as is implied in (fT3|) . and the 
errors cXi are model independent, what also is the case here. For the SNIa sample used 
we find x^{p) — X^iP^b/c) = ln(c/27r) ^ 7.2. 

To determine the best fit parameters for each model, we minimize x^{p), what is 
equivalent to maximizing the likelihood 

jC{p) oc e-^'(^)/l (18) 

The Bayesian evidence can then be calculated as 

E = j c{p)p{p)dp =Ypj^ (^^^ 

where P{p) is the prior probability distribution for the parameters, which we adopt to 
be flat, and Vp is the volume in the parameter space defined by the prior intervals. 
We chose conservative prior intervals (see table [1]) based on physical considerations and 
"prior" information, i. e. previous results obtained with older and smaller SNIa samples 
[22l [25| [27l [34] . The prior boundaries are chosen to be large enough so most of the 
likelihood is retained in the integration f[T^ . but not too large to do not excessively 
penalize the Bayesian evidence through Vp. For all the models considered here, except 
for M2, as will be discussed in the next section, the Sa boundaries in the likelihood are 
well inside the prior volume, so the Bayesian evidence decreases linearly with Vp. 

We are able to compare models by calculating the Bayes factor between any two 
models Mj and M^, which we define as 

B -^^^^^ (20) 

Some authors [35|, |36] offer qualitative interpretations of the Bayes factor value (Jeffreys 
scale) that say how one model is favoured over the other, given the data and priors. 
Note that under our convention of the Bayes factor, if E{Mj) > E{Mi) than In Bij is 
positive. 



4. Results and Discussion 



In figure [2] we show the likelihood results for the kinematic models considered and 
also for the ACDM model for comparison. The horizontal axis depicts the deceleration 
parameter today in the same scale for all models to facilitate the comparison among 
models. At the panels for ACDM and Mq the full likelihood is plotted as a function 
of go — we use that go = f^m — 1 for ACDM, see flA.22p — since these models have 
only one free parameter. For the other models it is presented the confidence contours in 
two dimensional parameter spaces (we marginalize over Zt for model M2 that has three 
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Figure 2. Likelihood results for all studied models (labels on the top of each panel). 
The top left panel is for models with one free parameter (ACDM and Mq). The 
remaining panels show the likelihood contours in the parameter space of each model, 
delimiting the la, 2a and 3a confidence regions. Model Mi is also shown in the (go Jo) 
parameter space. The dotted line in three of the panels depicts the image of ACDM 
models in the ((jo Jo) parameter space representation and the error bar corresponds to 
the la region for the best ACDM fit. The solid straight lines on the panels for ACDM 
and M3 depict the allowed regions in parameter space for these models. 
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-0.74 ±0.22 
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M4 : dL{ 


z) expansion 


-0.66 ±0.20 
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Table 1. Models, maximum likelihood estimates of the universal kinematic parameters 
(<Zo, Jo, -Zt) and Icr projected errors, prior volume, goodness of fit and Bayes factor in 
relation to Mq. All models are in a flat universe, Vp is the prior volume; priors are 
go £ [—2,0] for all kinematic models, qi G [—1,4] for Mi, qi G [—1,5] and zt G [0, 1.5] 
for M2, jo G [-5,1/8] for A/3, jo G [-5,1] for Nh, and a„ G [0,1] for ACDM. 
The relations between the particular model parameters and the universal kinematic 
parameters are given in the Appendix. 



degrees of freedom). The equations for M3 — see flA.15P in the Appendix — indicate 
that the maximum physical value for jo in this model is 1/8, so we limit the graph to 
this region. 

Table [1] contains the maximum likelihood values and la projected errors for the 
deceleration and jerk parameter today and for the transition redshift. Note that some 
of these parameters are the best fit values for free parameters of the models and others 
are derived quantities from them. See the Appendix for the expressions. Table [1] also 
presents the parameter space volume of the priors for each model; the goodness of fit 
as quantified by Xred = Xminl ~ ^p), where N is the number of data points and Up is 
the number of free parameters in each model; and the Bayes factor in relation to model 
Mo. 

The deceleration today is significantly negative in all models, but its exact value 
and uncertainty are mo del- dependent. The likelihoods for models with just one degree 
of freedom (Mo and ACDM) are very peaked in go (see top left panel of figure [2]), yielding 
a more precise determination of this parameter, and consequentially also of jo? than in 
the other models (see table [1]). For the models of constant acceleration today (Mo and 
M2) the current values of the cosmic jerk are in disagreement with the values found for 
the remaining models, which are compatible among themselves and with the ACDM 
prediction. The values of the transition redshift determined in all kinematic models are 
close to Zt ~ 0.5, and are compatible in lex with the higher value obtained in the ACDM 
model. 

Our likelihood contours for model Mi at figure [2] are qualitatively similar, but with 
tighter constraints, to what was obtained in previous works [3l [25]. We also represent 
the Ml model in the parameter space of the deceleration and jerk today, making use of 

(HI). 
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The panel for M2 shows a quahtatively similar plot to the last panel of figure 
2 of Shapiro & Turner [22]. The Zt marginalized likelihood contours for M2 at 2a 
and 3(7 confidence levels show that the data do not constrain strongly the deceleration 
parameter for redshifts above the transition. The parameter qi can assume very large 
positive values and also negative values. In this last case, represented by the points 
below the dashed line on the M2 panel of figure [21 2;^ is a transition redshift between two 
accelerated phases. Points above the dashed line represent situations in which there is 
a transition from a decelerated phase at high z to an accelerated phase at low z. We 
also tested a larger upper prior boundary for qi, but even for qi = 10 the 2a likelihood 
contour is open (in this case lni?02 = 1-6). 

Models M3 and M4 (bottom panels at figure [2]), and Mi as well, have roughly 
similar likelihood contours in the {qo,jo) plane due to the relation ([5]), implying similar 
determinations of go and jo, similar goodness of fit and Bayesian evidences. These 
likelihood contours also show that the best fit ACDM model, represented by a la error 
bar over the dotted line with jo = is compatible at la confidence level with Mi, M3 
and M4. Similar result was obtained by Rapetti et al. [27j using a joint analysis of the 
Gold [3] and SNLS [1] SNIa samples and X-ray cluster gas mass fraction measurements. 
When considering separately the SNIa sets, those authors found jo = —2.8til ^^e 
Gold set and jo = — l-3^i;4 for the SNLS set. We also examined the case where just 
supernovae with z < 1 (289 events) were fitted by model M4 and the results do not 
differ considerably from the ones for the full sample (307 events). In fact, the goodness 
of fit is worst for the subsample with z < 1, xled — 1-057. 

All models, except Mq, have similar goodness of fit, as quantified by Xred- fact, 
if we observe the best fit curves at figure [Tj it is very difficult to judge which one best 
describes the data. The Bayes factor segregates Mo from the other models, but only 
weakly disfavours M2 over the remaining models, mainly because of the penalizing effect 
of the larger prior volume of M2 (due to it having three free parameters) in relation to 
the other models. The basic result is that the Bayes factor is unable to significantly 
distinguish models Mi, M3, M4 and ACDM. This is a somewhat different result from 
what was found by Elgar0y & Multamaki [25] when examining separately the Gold [3] 
and SNLS ^ SNIa sets. There the authors obtained a clear ranking of the models based 
on the Bayes factor (see their table 5), even though these rankings were different for the 
two SNIa samples. 

It is worth to note that the behaviour for the evolution of the deceleration parameter 
is very distinct among the models considered here. The simplest kinematic models, 
Mq, Ml and M2, have self-evident q{z), but M3 has a more subtle, and interesting 
acceleration history. For this model q{z) — given by (]A.12p — has a qualitatively 
similar behaviour to that of ACDM (]A.22p . namely, asymptotic constant values in the 
past and future with a smooth transition around zt. For ACDM the asymptotic limits 
are fixed, q{z 00) = 0.5 and q{z —1) = —1, but for M3 they will depend on 
go and jo. In the specific case of our best fit values for M3, q{z 00) ~ 0.76 and 
q{z — > — 1) ~ —1.3. There is no deceleration in the future in this model, and therefore, 
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no new phase transition. For M4 the extrapolation of q{z) is not expected to be vahd 
far from the redshift range for which it was adjusted, and, in fact, flA.18|) presents 
nonphysical behaviour in the extremes past and future. None of our models are able 
to hint at any slowing down of the cosmic acceleration today or in the future, such as 
suggested recently [37] in the context of an extended SNIa sample. 

5. Summary and Conclusions 

Bellow we summarize our main accomplishments and conclusions. 

1. We perform, for the first time, a kinematic analysis of the 307 SNIa compiled 
in the Union set [31]. This approach, also called Friedmanless, allows us to analyse the 
cosmic expansion just having to assume the homogeneity and isotropy of the universe, 
and not having to make any assumption about the underlying gravitational theory and 
energy components of the universe. 

2. We employ several kinematic models that were used by different groups with 
various SNIa samples before. Using a unified framework and a single data set, we are 
able to compare the kinematic models directly. We calculate for each model the goodness 
of fit (as measured by xled) Bayes factor. Even very simplistic kinematic models 
can give an equivalent description of the cosmic expansion to the one provided by the 
currently favoured concordance ACDM model. More to the point, current data is not 
powerful enough to clearly discriminate among some of these simple models. This is 
a distinct conclusion from what Elgar0y & Multamaki [25] obtained using separately 
the SNLS and Gold samples for a particular class of the models studied here, and who 
found conflicting results from the two samples. Nevertheless, some kinematic models 
(Ml, M3 and M4) were shown to be superior, or more realistic, than kinematic models 
with constant deceleration (Mq and M2). 

3. We give, for all kinematic models studied, the expressions and estimates 
for a minimal set of parameters that characterizes the recent history of the cosmic 
expansion: go (deceleration today), Jq (cosmic jerk today) and Zt (transition redshift 
from a decelerated to an accelerated phase). 

4. Independently of models, the universe is in a phase of accelerated expansion, 
however the value of the deceleration parameter today is model-dependent. For the 
most realistic kinematic models (Mi, M3 and M4), go is in the la range [—0.96, —0.46]. 

5. There is evidence that the deceleration parameter was higher and positive in 
the past, implying a transition from a decelerated phase to an accelerated one. The 
transition redshift between these two phases is found to be around 0.5 in all kinematic 
models, being in the la range [0.36, 0.84] for the the most realistic kinematic models. 
That is compatible in la with the higher value predicted by ACDM. 

6. The value of the cosmic jerk today can be used as a measure of a possible 
deviation from a ACDM model, which exactly predicts jo = ~1 by definition. For the 
models with constant deceleration today (Mo and M2) the value for jo is significantly 
higher than for ACDM, but these models no longer have a strong cosmological appeal. 
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The remaining kinematic models (Mi, M3 and M4) have a shghtly preference for 
jo < —1, being in the la range [—3.2, —0.3]. But at the current confidence level yielded 
by the data there is no significant departure from the ACDM prediction. We note 
that the constant jerk model is an interesting parameterization of the cosmic expansion 
because, among other reasons, it contains popular models as ACDM and EdS. 

Finally, taking into account the discussion in the previous section, about the 
behaviour of the deceleration parameter in the context of jerk models, it could be 
interesting to extend the present work by allowing additional contributions from a snap 
parameter (dependent on the fourth order derivative of the scale factor). Either in a 
simplified model, s{z) = constant, as well as in a expansion of the luminosity distance. 
Hopefully, this may help us to have some indication about a possible dynamic transition 
in the future, with the universe entering in a new decelerating phase [38]. Some work 
along these lines will presented in a forthcoming communication. 
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Appendix A. Kinematic Expressions 

In this appendix we show the analytical expressions for the basic quantities in all models 
investigated in the present work. 



M, 







H{z) = H,{1 + z) 



l+qo 
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j{z) = -{qo + 2ql) 
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go, z < zt 
qi, z> Zt 
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From flA.lSP we see that jo < |- 

M4 - defined by the expanded luminosity distance ([71), dL{z) 
where ^4 = (1 — go)/2 and B = —(1 — go — 3go — jo)/6- 
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(A.9) 

(A.IO) 

(A.ll) 
(A.12) 
(A.13) 

(A.14) 

(A.15) 
(A.16) 

[z + Az^ + Bz"^), 

(A.17) 
(A.18) 



ACDM, Vtra 



l + 2Az + {A + 3B)z^ + 2Bz^ 
j{z) = -[q + 2q' + {l + z)q'] (A.19) 
Zt : the real root of 1-2A- 2{A + 3B)zt -{A + 9B)z^ - 2Bzf = 0(A.20) 

H{z) = Ho [nui + zf + (1 - n^)] ^ (A.21) 

q{z) = [(1 + zf - 2(l/a„ - 1)] / [2(1 + zf + 2(l/fi„ - 1)] (A.22) 

J{z) = -1 (A.23) 

Zt = ~ I (A.24) 

Note that the expressions for ACDM can be easily obtained from flA.imA.T4l) . putting 
jo = —1 into OA.lSp . If we do similarly for Mq, putting jo = —(go + 2go) into OA. 151) . 
we recover the right expressions for H{z) and q{z), and, curiously, Zt = —1. This just 
illustrates that all models with a particular constant jerk are particular cases of M3. 
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